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Abstract. The experimental data on hadron yields and ratios in central Pb+Pb 
and Au+Au collisions at SPS and RHIC energies, respectively, are analysed 
within a two-source statistical model of an ideal hadron gas. These two sources 
represent the expanding system of colliding heavy ions, where the hot central 
fireball is embedded in a larger but cooler fireball. The volume of the central 
source increases with rising bombarding energy. Results of the two-source model 
fit to RHIC experimental data at midrapidity coincide with the results of the one- 
source thermal model fit, indicating the formation of an extended fireball, which 
is three times larger than the corresponding core at SPS. 



1. Introduction. Two-source model 

Searching for the quark-gluon plasma (QGP) is one of the main goals in the study 
of relativistic heavy ion collisions. The principal question is whether the strongly 
interacting matter reaches the stage of chemical and thermal equilibrium. Although 
this idea was put forward by Fermi 50 years ago |j| , up to now there is no unambiguous 
test to probe the degree of equilibration in the system. One of the possible 
approaches is to study the equilibration process within microscopic models [^]-^| . The 
more traditional way is to fit yields and transverse spectra of particles obtained in 
experiments to the statistical model (SM) of a fully equilibrated hadron gas (see |^| |l4| 
and references therein). In the model the macroscopic characteristics of the system, 
such as the particle number density, the energy density, etc., are derived via a set of 
distribution functions (we assume c = h = ks = 1) 
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where p and are the momentum and the mass of the hadron species i, T is the 
temperature, \ib and f^s are t ne baryon chemical potential and strangeness chemical 
potential, Bi and Sj are the baryon charge and strange charge of hadron i. The sign 
"+" in Eq. dU) is for fermions and sign"—" for bosons. 

As the experimental data became more precise, it has been understood that 
the ideal SM does not provide an adequate description of all hadron multiplicities. 
Particularly, the yields of pions are usually underestimated while the abundances 
of strange particles are overpredicted. Therefore, some modifications to the SM 
have been proposed, such as excluded volume effects f§|, |9|; strangeness suppression, 
75 < 1 chemical non-equilibrium of light quarks The assumption of a 

single expanding source is the basic ad hoc hypothesis of these models. However, if 
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the baryon density and/or the strangeness density are not the same everywhere in the 
reaction volume, the scenario with several independent sources cannot be reduced to 
the single source scenario. Our investigations of the two-source scenario are inspired 
by the experimental observation fTq ] and microscopic model predictions Q of low net 
baryon densities in the midrapidity range of relativistic heavy-ion collisions at SPS 
energies. Thus, local equilibrium may occur separately in different zones of the fireball. 

The proposed two-source statistical model (TSM) of a hadron gas |l6|] divides the 
whole reaction zone into two regions: the outer region (source 1 or SI) and the inner 
region (source 2 or S2). The two sources are allowed to possess different temperatures, 
net baryon and strangeness densities, etc. The characteristics of each of the fireballs 
can be described via the four independent parameters, such as volume V, fireball 
temperature T, net baryon density ps, and net strangeness density ps- Compared to 
the SM, the number of free parameters in the two-source model increases to seven: 
Although the net strangeness in each of the sources can be nonzero, they are linked 
via the the condition of total strangeness conservation N$i + N$2 = , i.e., either psi 
or ,052 may be considered as a free parameter. 

2. Hadron production at SPS and RHIC 

The baryon yield and ratios of hadrons at midrapidity in central Pb+Pb collisions at 
158 AGeV are listed in Table together with the results of the TSM and the SM fit. 
All hadrons with masses less than 2 GeV/c 2 are included in the fitting procedure. No 
additional constraints such as strangeness suppression or excluded volume are assumed 
except the feeding-back effect from resonance decay. Compared to the ideal SM, 
the TSM improves the agreement with the experimental data. The thermodynamic 
quantities obtained from the two fits to experimental data are shown in Table |l| also. 
One can see that the two-source object can be interpreted as a hot, relatively small 
core surrounded by a cooler and larger halo. The major part of baryons is contained 
in the outer source, while the inner source contains almost all antibaryons. 

The strangeness density is negative in S2 and positive in SI, i.e., the inner source 
contains more s-quarks than s-quarks. This finding is supported by microscopic model 
calculations [Q. The possible explanation of the phenomenon is as follows: strange 
and anti-strange particles must be produced in pairs. Because of the small interaction 
cross section with hadrons, K + and K° are leaving the central reaction zone easier 
than strange particles which carry s quarks, A and K 1 thus leading to a negative 
strangeness density in the midrapidity range. 

The energy density in S2 is about three times larger than that in SI. Such a 
low energy density in the outer source corresponds to the energy density at thermal 
freeze-out rather than at chemical freeze-out 0. In other words, the solution for two 
sources cannot be reduced to the one-source picture even in the case where exclusively 
midrapidity data have been used. Similar fit has been performed to solely 47r-data 
and to a mixture of 47r-data with midrapidity data. Results of all three fits favour the 
idea of the formation of a compact hot baryon-dilute central zone with the following 
averaged characteristics: temperature T = 157 ± 2 MeV, volume V = 0.6 ± 0.1 Vq, 
where Vq is the volume of a lead nucleus, and baryon chemical potential ps = 31 ± 14 
MeV. The temperature of the halo is much lower, Ts 1 = 117 ±3 MeV. 

Experimental data on hadron yields and ratios in the midrapidity range in central 
Au+Au collisions at y^i = 130 AGeV are listed in Table || together with the predictions 
of the SM and TSM. Surprisingly, now the results of the SM fit and the TSM fit are 
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Table 1. Baryon yield and hadron ratios at midrapidity for central lead-lead 
collisions at SPS energies and predictions of the single-source and two-source 
statistical models of an ideal hadron gas. 





Data 


SM 


TSM 


Coll. 




SI 


S2 


SM 


N B {net) 


372±10 


371.9 


372.1 


NA49 


T[MeV] 


117 


155 


158 


K+/K- 


1.85± 0.1 


1.97 


1.87 


NA44 


V[fm 3 ] 


7250 


1705 


4203 


P/P 


0.07± 0.01 


0.069 


0.058 


NA44 


N B 


348 


55 


408 




0.249± 0.019 


0.231 


0.247 


WA97 


N b 





30 


36 


n/n 


0.383± 0.081 


0.427 


0.405 


WA97 


N s 


34 


42 


117 


A/A 


0.128± 0.012 


0.130 


0.137 


WA97 


N s 


21 


55 


117 


n/n 


0.081±0.013 


0.133 


0.108 


WA98 


V B [MeV] 


460 


45 


213 




0.125±0.019 


0.121 


0.120 


NA49 










K° s /h- 


0.123±0.02 


0.102 


0.102 


WA97 










A/h- 


0.077±0.011 


0.069 


0.064 


WA97 










Q/E 


0.219±0.045 


0.131 


0.104 


WA97 










S-/A 
X 2 /DOF 


0.110±0.01 


0.156 


0.115 


WA97 












46/9 


16/5 













almost identical. It seems that the volume of the central fireball significantly increases 
by the transition from the SPS energies to the RHIC ones, and that hadrons detected 
in the midrapidity region are originated from a single thermalized source. Its volume 
is more than 5000 fm 3 , i.e., three times larger than the core volume at SPS, and 
the temperature reaches 186 MeV. However, if the multiplicity of negatively charged 
hadrons, h~ , is excluded from the data set, the temperature of the modelled system 
drops to 176 MeV, that is extremely close to the value T — 175 MeV obtained in |Q. 
This important question should be clarified in future studies. We checked that the 
incorporation of the excluded volume effects by assigning the hard-core radius r = 0.4 
fm to all particles leads to an enlargement of the volume but does not affect the 
temperature of the fireball. 

Another useful characteristics at the chemical freeze-out are the energy per 
particle and the entropy per baryon. In jl0| the criterion E/N s» 1 GeV was 
introduced. The predictions of the SM and the TSM for the midrapidity range are as 
follows: E/N = 1.1 GeV, and S/A = s/p B = 121. The last value is about 20% below 
the value s/ps = 150 predicted by the UrQMD calculations, but the latter are for 
Au+Au collisions at y/s — 200 AGeV J!?]]. It would be interesting to perform these 
microscopic model calculations also at y^i = 130 AGeV. 

3. Conclusions 

In summary, the TSM fit to the experimental data taken at midrapidity at RHIC 
coincides with the standard single-source fit. This result supports the idea of a 
formation of an extended hot fireball in the central zone of heavy-ion collisions at RHIC 
energies. The temperature of the central fireball varies from 176 MeV to 185 MeV, 
depending on the incorporation of the multiplicity of negatively charged hadrons in 
the fit, whereas the excluded volume effects seem not to affect the fireball temperature. 

At SPS it is found that the properties of the system at chemical freeze-out can 
be well understood in terms of two sources, a central core and a surrounding halo, 
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Table 2. 


The same 


as Table 


but for Au+Au collisions at ^/s = 


130 AGeV. 






Data 


SM 


TSM 


Coll. 




SM 


SI 


S2 


Nb (net) 


343±11 


340.6 


340.6 


PHOBOS 


T[MeV] 


185.7 


185.9 


185.5 


hr 


2050±250 


2238 


2239 


PHOBOS 


V\fm 3 ] 

V [J III J 


5157 


2608 


2560 


p/p 


0.60± 0.07 


0.57 


0.57 


PHOBOS 


N B 


893 


456 


437 


p/iT- 


0.08± 0.01 


0.087 


0.087 


STAR 


% 


553 


281 


272 


K-/K+ 


0.87± 0.08 


0.75 


0.75 


STAR 


N s 


307 


155 


152 


K-/-K- 


0.149±0.02 


0.153 


0.153 


STAR 


Ns 


307 


155 


152 


K*°/h- 


0.060±0.017 


0.036 


0.036 


STAR 


fiB [MeV] 


53 


53 


52 


K*°/h~ 


0.058±0.017 


0.030 


0.030 


STAR 










A/A 


0.77± 0.07 


0.69 


0.69 


STAR 












0.82± 0.08 


0.80 


0.80 


STAR 










X 2 /DOF 




9.6/3 


9.6/7 













both being in local chemical and thermal equilibrium. Temperatures as well as baryon 
charge and strange charge of the two sources are different. 

Strangeness seems to be in equilibrium in both sources. This observation is in 
line with the fact that js ~ 1 if one fits the particle ratios from the midrapidity range 
to the SM ||, but 75 < 1 if one intends to fit 47r-data A possible explanation 

for this puzzle is a non-homogeneous distribution of the strange charge within the 
reaction volume. It would be interesting to study the forthcoming 4-7r-data on Au+Au 
collisions at RHIC energies in order to check (i) the increase of the volume of the 
central fireball w.r.t. the halo; (ii) equilibration of strangeness in both sources; (iii) a 
possible change of the halo temperature. 

Acknowledgements. We are grateful to L. Bravina for fruitful discussions. The 
work was supported by the DFG, BMBF under the contract No. 06TU986, and the 
NSF of China under the contracts No. 19975075 and No. 19775068. 

References 

[1] Fermi E 1950 Prog. Theor. Phys. 5 570 

[2] Geiger K and Kapusta J I 1993 Phys. Rev. D 47 4905 

[3] Bravina L V et al 1998 Phys. Lett. B 434 379 

Bravina L V et al 1999 J. Phys. G: Nucl. Phys. 25 351 

[4] Bravina L V et al 1999 Phys. Rev. C 60 024904 

Bravina L V et al 2000 Phys. Rev. C 62 064906 

[5] Sollfrank J, Heinz U, Sorge H and Xu N 1999 Phys. Rev. C 59 1637 

[6] Rafelski J 1991 Phys. Lett. B 62 333 

[7] Becattini F 1996 Z. Phys. C 69 485 

[8] Yen G D and Gorenstein M I 1999 Phys. Rev. C 59 2788 

[9] Braun-Munzinger P, Heppe I and Stachel J 1999 Phys. Lett. B 465 15 

[10] Cleymans J and Redlich K 1999 Phys. Rev. C 60 054908 

[11] Letessier J and Rafelski J 1999 Phys. Rev. C 59 947 

[12] Lu Z D et al 1998 High Energy Phys. Nucl. Phys. 22 910 

[13] Zhang Z Y et al 2000 Phys. Rev. C 61 065204 

[14] Braun-Munzinger P, Magestro D, Redlich K and Stachel J 2001 Phys. Lett. B bf 518 41 

[15] Ambrosini G et al 1998 Phys. Lett. B 417 201 

[16] Lu Z D, Faessler A, Fuchs C and Zabrodin E 2001 Preprint friucl-th/011004C 

[17] Bravina L V et al 2001 Phys. Rev. C 63 064902 

Bravina L V et al 2001 J. Phys. G: Nucl. Phys. 27 421 



